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PULSE  COMPRESSION  DEGRADATION  DUE  TO  OPEN  LOOP 
ADAPTIVE  CANCELLATION,  PART  II 


1.  INTRODUCTION 

An  exact  expression  for  the  perturbed  sidelobe  level  of  a  compressed  pulse  that  has  been  pre- 
processed  through  an  adaptive  canceller  was  derived  in  Ref.  1.  The  pertinent  assumptions  of  that 
analysis  are: 

1.  the  adaptive  canceller  was  implemented  by  using  the  Sampled  Matrix  Inversion  (SMI)  algo¬ 
rithm  [2]  or  its  equivalent,  the  Gram-Schmidt  canceller  [3], 

2.  the  input  noises  were  temporally  independent  and  Gaussian, 

3.  the  desired  signal’s  input  vector  (or  code)  was  completely  contained  within  the  samples  that 
were  used  to  calculate  the  adaptive  weights  and  is  only  present  in  the  main  channel,  and 

4.  the  adaptive  weights  were  computed  from  the  same  data  set  to  which  they  are  applied  (con¬ 
current  processing). 

Earlier  research  has  shown  that  because  of  finite  sampling,  the  quiescent  compressed  pulse 
sidelobe  levels  are  degraded  by  preprocessing  the  main  channel  input  data  stream  (the  uncompressed 
pulse)  through  the  adaptive  canceller.  It  was  also  shown  that  the  level  of  degradation  is  independent 
of  whether  pulse  compression  occurs  before  or  after  the  adaptive  canceller  under  assumption  three. 

The  exact  expression  [1]  for  pulse  compression  degradation  requires  computer  assistance  to 
evaluate  this  expression.  This  report  derives  a  “rule  of  thumb”  expression  that  is  a  good  approxima¬ 
tion  to  the  exact  expression. 

2.  BACKGROUND 

Figure  1  is  a  functional  block  diagram  of  an  adaptive  canceller  followed  by  a  pulse  compressor. 
The  adaptive  canceller  linearly  weights  the  auxiliary  channels  with  weights  that  are  calculated  from  a 
batch  of  sampled  input  data.  The  main  channel  consists  of  a  desired  signal  plus  noise  that  may  or 
may  not  be  correlated  with  the  auxiliary  channels.  It  was  shown  in  Ref.  1  that  when  analysing  the 
pulse  compression  degradation  it  is  necessary  only  to  consider  the  interaction  of  the  main  channel's 
desired  signal  with  the  random  variables  in  the  auxiliary  channels  (Fig.  1).  Thus  for  analysis  pur¬ 
poses,  the  adaptive  weights  of*,,,  n  =  1,  2,  ...,  N  -  1  are  only  a  function  of  the  desired  signal  s 
and  the  samples  of*,,.  Furthermore,  as  the  number  of  independent  samples  goes  to  infinity,  the  auxi¬ 
liary  adaptive  weights  go  to  zero  [1]. 


Manuscript  approved  April  4.  1991. 


KARL  G12RLACH 


Fig.  1  —  GS  canceller  followed  by 
a  matched  Filter 


In  Fig.  1,  s  represents  the  desired  signal  vector  (or  code)  of  length  L,  and 
x„,  n  =  1,2,  ...,  N  -  1  represents  the  Mh  auxiliary  random  data  vector  of  length  K.  The  canceller 
shown  is  the  Gram-S~hmidl  (GS),  which  is  numerically  equivalent  to  the  SMI  algorithm  [3].  We 
denote  it  by  GSK  ^  where  K  is  the  number  of  samples  per  channel  used  to  calculate  the  canceller 
weights  and  N  is  the  number  of  input  channels  (main  and  auxiliaries). 

The  pulse  compressor  is  essentially  the  matched  Filler  for  a  given  radar  waveform.  Most  of  the 
energy  in  the  received  radar  waveform  is  compressed  into  a  given  single  range  cell  and,  thus,  the  sig¬ 
nal  level  can  be  increased  significantly  for  detection  purposes.  However,  some  energy  does  leak  into 
the  sidelobes  of  the  compressed  pulse  response,  resulting  in  low  gain  in  range  cells  outside  of  the 
given  range  cell.  If  a  target  or  piece  of  clutter  is  large  enough,  it  can  break  through  and  be  detected 
in  these  range  sidelobes,  falsely  indicating  a  target  detection  or  masking  a  real  target.  Thus  it  is 
highly  desirable  to  maintains  low  sidelobe  response. 

Let  r  equal  the  2 L  -  1  output  vector  of  the  pulse  compressor.  If  an  adaptive  canceller  is  not 
being  used,  then  it  is  straightforward  to  show  that 

r  =  S's  (1) 
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where 


s  —  (5i.  i-2>  »  s,)r, 
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and  T,t  denotes  transpose  and  complex  conjugate  transpose,  respectively.  S  is  a  L  x  {2L  -  1) 
matrix  called  the  autocorrelation  function  (ACF)  matrix  of  s.  If  L  <  K,  we  define  an  augmented  sig 
nal  vector  s„„A,  of  length  K  such  that  the  first  L  elements  are  elements  of  s  and  the  remaining  elements 
are  zero.  Sm,s  is  defined  as  the  augmented  (2 K  -  1)  x  (2 K  -  1)  ACF  matrix  of  s  using  the  ele¬ 
ments  of  sailg.  The  quantity  rallg  is  defined  as  the  augmented  2 K  -  1  output  vector  of  the  pulse 
compressor.  Thus 


«“8 


-•  C'  c 
Uaug°migm 


(3) 


Let  s'  be  the  resultant  output  vector  after  s  has  been  processed  through  the  GS  canceller  and  s,;„A>  be 
the  resultant  augmented  GS  output  vector.  This  resultant  output  vector  is  then  inputted  to  the  matched 
filter  of  the  vector  s,  or  equivalently,  sailg.  If  we  set  r'  equal  to  the  response  of  s„'„A,  match  filtered 
with  saug  then 


r'  = 

*  'Jaug*augm 


(4) 


In  Ref.  1  it  was  shown  under  assumptions  I  through  4  (given  in  Section  1)  that  the  average 
pulse  compressed  sidelobe  level  after  adaptive  cancellation  is  given  by 


SLa(l)  = 


(K  -  N  +  1  )(K  -  AM-  2) 


SLJl)  + 


K(K  +  1) 


(K  -  N  +  \)(K  -  N  +  2) 


/1i2(^.A')||sc(/)||2,(5) 
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where 

SLa(I)  is  the  average  pulse  compressed  bidelobe  level  after  adaptive  cancellation  of  the  /th  range 
sidelobe  (sidelobes  are  numbered  ±1,  l  =  1,  2,  these  can  be  related  directly  to  the 
elements  of  r';  for  example,  /  =  ±  1  are  the  sidelobes  adjacent  to  the  match  point). 

SLp  l)  is  the  quiescent  pulse  compresbed  sidelobe  level  of  the  /th  sidelobe  (K  =  oo  or  equivalently 
no  adaptive  cancellation  before  pulse  compression,  these  can  be  related  directly  to  the  ele¬ 
ments  of  r) 

K  is  the  number  of  independent  samples  per  channel  used  to  calculate  the  adaptive  canceller 

weights 

N  is  the  number  of  channel  (main  and  auxiliaries) 

sc(/)  is  the  K-hh  column  of  the  augment  AFC  matrix  Snug,  l  =£  K,  and 

IM0!i?  =s'(/)sc(/) 

Note  that  SL„(l)  and  SLtl(l)  are  normalized  to  the  mainlobe  pulse  compression  gain  (adapted  or  quies¬ 
cent.  respectively)  which  is  set  equal  to  one  or  0  dB. 

The  scalars  AU(K,N)  and  A  12(AT, Af)  are  computed  as  follows.  Consider  the  two  parallel  adap¬ 
tive  cancellers  shown  in  Fig.  2.  Define 

u0,v0  are  arbitrary  AT-length  main  channel  input  vectors, 

are  AMength  main  channel  output  vector,  and 

x„  =  (.r„(i),  .v„(2),  ...,  x„(K)f,  n  =  1,  2,  ...,  N  -  1,  ^-length  random  data  vector 
of  the  nth  auxiliary  channel. 

The  elements  of  x„,  n  =  1,2 . /V  --  1,  are  assumed  to  have  the  following  characteristics: 

1.  x„(k),  n  =  1,  ...,  N  -  \,  k  =  \,  ....  K  are  identically  distributed  circular  Gaussian  com¬ 
plex  random  variables  (r.v.) 

2.  E\xn(k)\  =  0,  E\\x„(k,  j2|  =  1,  where  £{  |  denotes  expectation  and  |  |  denotes  magni¬ 
tude 

3.  E\x„i(ki)xt(k2)\  =  0  unless  =  n2  and  k |  -  k2,  where  *  denotes  complex  conjugate. 
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Fig.  2  —  Parallel  AMnput  GS  cancellers 
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^{K+iflK+iU2}  =  £{Kv„|2} 


_ 1 _ 

(K  -  n)(K  -  n  +  1) 


+  £(K||2!M|2} 


2  ,  1 _ 

K  —  it  '  (K  -  n)(K  -  n  +  1) 


(ID 


where//  =  0,  1.  N  —  1. 

3.  BOUND  DERIVATION 

The  expression  derived  for  SL„  given  by  Eq.  (5),  although  exact,  does  not  readily  indicate  how 
the  adaptive  sidelobe  level  varies  with  N  and  K.  In  this  section  we  derive  a  tight  upper  bound  on  SLa 
that  is  in  terms  of  explicit  expressions  of  K  and  /V. 

This  bound  is  obtained  by  considering  Eqs.  (10)  and  (11).  Instead  of  deriving  a  recursive  rela¬ 
tionship  for  £(SJu„!|2j|v„][2}  in  this  equation,  we  upper-bound  this  expectation  by  using  the  inequality 
[4] 


£fl|uj2||v,,||2)  ^  VillM4}  £{||v„|j4}.  (12) 

This  inequality  is  merely  another  form  of  the  Cauchy-Schwartz  inequality.  It  allows  us  to  upper- 
bound  the  joint  moment  in  terms  of  moments  of  individual  random  variables. 

It  was  shown  in  Ref.  3  that 

£||Kli4)  -  IK  II4  iK  ~  "*f  (13) 

and 

£!l|v„llJi  =  l|v„||J  <M> 

Substituting  these  expressions  into  Eq.  (12)  results  in 

cilKlPlWpi  s  IK  ||2  IK  |p  fg  *  '*•  (l5) 


Substituting  Eq.  (15)  into  Eq,  (10)  results  in 


£[K,-|V„+|  |2|  s  E\  |u'v„ 


_ I _ 

( K  -  n)(K  -  n  +  1) 


+  II»o!!2!|v0||2 


I 

K{K  +  1)' 


(16) 
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It  is  apparent  from  Eq.  (16)  that  E[  |  u'v„  1 2),  n  —  0,  1,  ...,  N- 1  can  be  upper-bounded  by 
oi„  where  is  found  by  the  recursive  relationship 


W„4.|  =  1  - 


K  -  n  (K  -  n){K  -  n  +  1) 


Kir-iivoH2 

K(K  +  1) 


Initial  condition  (IC)  w0  =  I  uovo  1 2- 


We  can  show  that  coA'-i  has  the  form 


N-2  2 

wa'-i  =  II  1  ~  + 

»i=0  K  11 


K  -  n  +  (K  -  n)(K  -n  +  1)  'U°V°' 


+  c  HuolPHvoil2,  (18) 

where  c  is  a  constant  to  be  determined.  In  fact,  c  does  notdepend  on  the  initial  condition,  a  fact  that 
we  use  to  find  c.  For  u0  =  v0,  it  follows  from  Eq.  (13)  that 


<■>*-1  =  =  I 


Substituting  Eq.  (19)  into  Eq.  (18)  and  solving  for  c, 
_  2(AT  -  1)  .  (N  -  l)N  N~2 


K  K{K  +  1) 

It  is  shown  in  the  Appendix  that 


K  —  n  (a  —  n)(K  —  n  +  1) 


K  —  n  (a  -  n)(K  -  //  +  1) 


(K  -  at  +  2)07  -J) 

K(K  -  \)(K  +  I)  ' 

Thus  inserting  these  inequalities  into  Eq.  (19)  results  in 
**  ~\  2 

mil  (•>,  „  ,  N  —  \  ,  ,  ti  .  (E  —  N  +  2 )(<V  —  I)  „  Mi ,, 

£(kv-,v,v_,|  j<  i - -j7  I  »ov0 !  *  +  +-jj  ■  IKIMIvoih  (23) 

J 

To  find  a  bound  on  the  adaptive  compressed  sidclobe  level  as  was  done  in  Ref.  I.  we  set 
”0  =  sIIWR  and  v0  =  sf  where  s„„x  is  the  augmented  /if- length  signal  vector  and  s,  is  a  column  of  the 
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augmented  signal  matrix.  Wc  note  that  ||s||2  =  1  and  that  SLq  =  |s's|2.  It  was  shown  in  Ref.  1 
that  the  expected  value  of  the  match  point  of  the  compressed  pulse  preproccssed  b>  a  CSK  N  canceller 
is 


£(|s"s|2} 


(K  -  N  +  l)(K  -  N  ■ f  2) 
K(K  +  1) 


(24) 


Thus  if  we  divide  both  sides  of  Eq.  (23)  by  the  expected  value  of  the  match  point  we  find 


sljd  < 


i  - 


N  -  1 


K(K  -  N  +  2) 


SLJl)  +  [|sc 


i>2 


N  -  1 


(K  -  N  +  l)(K  -  1)  ‘ 


(25) 


We  set 


Q(K,N)  =  1  - 


N  -  1 


K(K  -  N  +  2)  ’ 


and 


Thus 


ASL,AK,N)  = 


N  -  1 


(K  -  N  +  1  )(K  -  1) 


SLJl)  <  Q(K.N)SLJl)  +  ASLJK,N)\\sc(l)\\2. 

Similarly,  define  the  quiescent  sidclobe  level  factor 

K(K  +  l)An(K,N) 


Q(K,N)  = 


(K  -  N  +  1  )(K  -  N  +  2) 


and  the  adaptive  side  perturbation 


ASLJK.N)  = 


K(K  +  1  )A  l2(K,N) 


{K  -  N  +  1  )(K  -  N  +  2) 

so  that  Eq.  (5)  can  be  rewritten  as 

SLJl)  =  Q(K.N)SLJl)  +  ASLJK.N)  [M/)I|2. 


(26) 

(27) 

(28) 

(29) 

(30) 

(31) 


4.  RESULTS 

Wc  now  demonstrate  in  graphical  form  that  Q(K,.X)  and  ASLJK. .Vj  arc  close  approximations  of 
the  quiescent  sidclobe  level  factor  Q(K.X)  and  the  adaptive  sidclobe  perturbation  ASLJK.X).  respec¬ 
tively.  Define  the  following  ratios 


Q(K.N) 
QlK.N) ' 


(32) 
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and 

ASLa(KrN) 

r  a  =  — = - .  (33) 

&Sla(K.N) 

Wc  set  Nattx  =  N  -  1  and  K  =  MN  where  M  is  a  positive  integer  and  calculate  rq  and  rs  vs 
Naax  and  M.  Wc  restrict  M  >  2.  Many  cases  were  run  (M  <  10.  N  <  100).  and  the  two  ratios 
were  always  less  than  one  and  lower-bounded  by  the  case  when  M  =  2.  Thus,  we  only  present  the 
curves  for  M  -  2.  The  close  approximation  is  verified  by  the  plots  of  rq  and  r ±  shown  in  Figs.  3 
and  4.  respectively.  The  worst-ease  approximation  of  Q{K.N )  by  Q(K.N)  occurs  when  :\aIU  ~  1. 
M  =  2.  In  this  case  rq  (dB)  =  — 1.76  dB. 


Fig.  3  —  ry  vs  M  ■  2 


An  even  better  approximation  of  Q[K.N)  was  found  by  using  the  expression 


ASLAK.N)  = 


iV  -  1 

CAT  -  iV  -F  1) K' 


Note  that  the  difference  between  the  expression  for  ASI,tl  given  by  Hqs.  (34)  and  (27)  is  that  the 
K  -  l  is  replaced  by  K.  Define  the  ratio 


ASUJI) 

ASLJI)' 


(35) 


Figure  5  plots  this  ratio  for  M  -  2  vs  A',nif .  Note  that  (he  worst-ease  approximation  occurs  when 
iV(BH  =f  2.  M  =  2.  In  this  case.  r'A  (dB)  =  -.51  dB. 
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Because  Q(K,  N)  and  A SLa(K,N)  are  close  approximations  to  Q(K,  N)  and  ASLa(K,N)  respec¬ 
tively,  it  is  straightforward  to  show  that 


An(K,N)=  1- 


N  -  1 


An(K,N) 


±  (K  -  N  +  2 )(N  -  1) 
K(K  -  1)(£  +  1) 


Again,  if  we  replace  K  —  1  with  K  in  Eq.  (37),  an  even  better  approximation  results: 

Aa(K,N)  =  (K  +  2){N  ~  U  .  (38) 

K2(K  +1) 

Inserting  the  approximate  expressions  given  by  Eq.  (26)  for  Q(K,N)  and  Eq.  (34)  for 
ASLa(K,N)  into  Eq.  (31)  results  in 


1  - Tur^h 2)  +  g H si  w 


Define  K2M(l)  to  be  the  minimum  number  of  independent  samples  such  that  SLa(l)  <  2 SL(I,  where 
SLq  =  max  SLq(l );  (i.e.,  the  average  adaptive  sidelobe  level  at  a  specific  range  sidelobe  /  is  at  most  3 

dB  above  the  maximum  quiescent  sidelobe  level).  It  is  straightforward  to  show  that 


^3dB(0  -  —  + 


N  -  1  “  IMOII2 

-  +  {N  -  1)  " - 

2  J  2 SLq  -  SLq(l) 


when  SLg(l)  «  1.  The  actual  number  of  samples  used  to  ensure  that  all  adaptive  sidelobes  are 


below  2 SLq  would  be 


^3dB  -  roax  A^bCO- 
U*  o 


If  the  maximum  quiescent  sidelobe  level  occurs  close  to  the  main  lobe,  then  j[s<  |p 
that  AT3di)(/)  is  maximized  at  this  maximum  quiescent  sidelobe  level.  Hence, 


=  1  and  we  find 


K  -  N  ~  1  + 
A3clB  -  2 -  + 


N  -  1 
2 


+  N^± 


Reference  1  pointed  out  that  the  pulse  compression  degradation  analysis  can  be  applied  to  quan¬ 
tifying  the  canceller  degradation  caused  by  a  desired  signal’s  presence  in  the  samples  used  to  calculate 
the  adaptive  weights.  If  the  desired  signal  has  the  power  o ]  after  pulse  compression,  then  the  average 
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power  residue  caused  by  the  signal  in  the  K  -  1  range  bins  not  containing  the  signal  can  be  shown  to 
equal  at  most  o2sASLa(K,N)  plus  possibly  the  signal  power  due  to  the  quiescent  compressed  sidelobes. 
Let  a2,in  be  the  quiescent  output  noise  power  level  of  the  canceller.  Define 

2 

8(1)  =  A SLa(K,N)  ||sc(/)||2  (43) 

°min- 

and 


8  =  rnax  8(1). 
/./*  o 


(44) 


One  normally  desires  8  <  1,  otherwise  the  desired  signal  generates  more  range  sidelobe  power  than 

the  noise  power  residue.  Because  max  ||sc||2  =  1  and  using  the  good  approximation  for  ASLa(K,N) 

i.i 

given  by  Eq.  (34)  then 


8  = 


N  -  1 

(K  -  N  +  l)K 


ffmin 


(45) 


It  is  desirable  to  know  the  number  of  independent  input  samples  A^0  such  that  8 
shown  that 


It  can  be 


(46) 


We  note  that  a2 /a2 „„  equals  the  output  signal-to-noise  power  ratio  (S/N)0llt  of  the  adaptive  canceller. 
Thus  Eq.  (46)  reduces  to 


Ko 


+ 


+  (N  -  1) 


Ollt 


(47) 


5.  SUMMARY 

An  exact  expression  for  the  perturbed  sidelobe  level  of  a  compressed  pulse  that  has  been  pre- 
processed  through  an  adaptive  canceller  was  derived  in  Ref.  1.  The  exact  expression  requires  com¬ 
puter  assistance  to  evaluate  this  expression.  In  this  report,  a  “rule  of  thumb’’  expression  is  derived 
that  is  a  good  approximation  to  the  exact  expression.  Furthermore,  this  same  approximation  can  be 
used  to  derive  a  good  approximation  for  the  canceller  noise  power  level  that  is  induced  by  having  a 
desired  signal  present  in  the  canceller  weight  calculation.  An  expression  for  the  number  of  indepen¬ 
dent  samples  necessary  to  equalize  the  signal-induced  power  with  the  quiescent  interference  level  is 
also  derived. 
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Appendix 

PROOF  OF  EQS.  (21)  AND  (22) 

We  first  prove  Eq.  (21).  Define 


b  =  TT  1 - ±—  +  - i _ 

„=o  I  K  —  n  (K  -  n)(K  -  n  +  1) 


—  |  __  ~  _ _ Z _  J 

11  =o  [_  K-Ii  (K  -  n)(K  -  n  +  1)  1  (K  -  n)(K  -  n  -  1) 


N-2 

-  n  i 

H=0 


2  _ 2 _  ^-2 

*  -  n  (K  -  n)(K  -  n  +  1)1  1 


(K  -  n)(K  -  n  -  I) 


n  1- 


K  -  n  (K  -  n)(K  -  n  +  1) 


=  (K  -  N  +  1  )(K  -  N  +  2) 
K{K  +  1) 


l  _  2(A  -  1)  A(/V  -  1) 
£  w  +  i) ' 


b  =  (K  —  N  +  1)(AT  -  N  +  2)  ""2  \ 

*c*  +  »  5  1 


(*  -  n){K  -  n  -  1) 


<  (K  -  N  +  l)(K  -  N  +  2)  "  "  2  I  1 

*(*  +  0  /A 


(K  ~  n)2 


^tt2  ] _ 1  _  (K  +  1)(AT  -  A  +  I) 

n  =o  [  (K  -  «)2J  K(K  -  N  +  2) 


b  <  ~  M  +  l)2  = 

a:2 


n  -  i 
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Next  we  prove  Eq.  (22).  Set 


Now 


N-2 


a  =  n 

n  =0 


( K  -  n)(K  -  n  -  1) 


and 


b  = 


,  +  2(iV  ~  1)  N  (N  -  1) 

K  K[K  +  1) 


a 


or 


a  > 


N-2 


n 

«=0 


1 

(K  -n  -  l)2' 


K(K  -  N) 

(K  -  1  )(K  -  N  +  1) 


!—«<!  — 


K(K  -  N) 

(K  -  l)(K  -  N  +  1) 


N  -  1 

(K  -  l)(ff  -  N  +  1)  ‘ 


Using  Eq.  (A7)  and  Eq.  (20),  it  can  be  shown  that 


c 


\  _  2  (N  -  1) 
K 


,  (A  ~ 

K(K  +  1) 


(1  -  a). 


(A6) 


(A7) 


(A8) 


(A9) 


(A10) 


Using  Eq.  (A9)  it  follows  that 


c  < 


1 


2  (N  -  1)  (/V  -  l)Af 

/T(AT  +  ,1) 


W  -  1 

(AT  -  !)(£  -  N  +  i) 


(All) 


or 


{K  ~  V  +  2)(Af  ~  1) 
OT  -  l)(tf  +  1) 


c  < 


(A  12) 


